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ABSTRACT 

We present the results of an interferometric study of 38 millimeter-wave lines of 
12 CH30H in the vicinity of the massive star forming region W3(OH/H20). These lines 
cover a wide range of excitation energies and line strengths, allowing for a detailed 
study of excitation mechanisms and opacities. In this paper we concentrate on the 
region around the water maser source W3(H20) and a region extending about 30 arcsec 
to the south and west of the hydroxyl maser source W3(OH). 

The methanol emitting region around W3(H20) has an extent of 2.0 x 1.2 arcsec 
(4400 x 2600 AU). The density is of order 10 7 cm -3 , sufficient to thermalize most of 
the methanol lines. The kinetic temperature is approximately 140 K and the methanol 
fractional abundance greater than 10 -6 , indicative of a high degree of grain mantle 
evaporation. The W3(H20) source contains sub-structure, with peaks corresponding 
to the TW source and Wyrowski's B/C, separated by 2500 AU in projection. The 
kinematics are consistent with these being distinct protostellar cores in a wide binary 
orbit and a dynamical mass for the region of a few tens of Mq. 

The extended methanol emission to the southwest of W3(OH) is seen strongly only 
from the lowest excitation lines and from lines known elsewhere to be class I methanol 
masers, namely the 84.5 GHz 5_i-4o E and 95.2 GHz 8q— 7\ A + lines. This suggests that 
this region, like class I maser sources, is dominated by collisional excitation. Within 
this region there are two compact clumps, which we denote as swA and swB, each 
about 15 arcsec (0.16 pc projected distance) away from W3(OH). Excitation analysis 
of these clumps indicates the presence of lines with inverted populations but only weak 
amplification. The sources swA and swB appear to have kinetic temperatures of order 
50-100 K and densities of order 10 5 -10 6 cm -3 . The methanol fractional abundance for 
the warmer clump is of order 10~ 7 , suggestive of partial grain mantle evaporation. The 
clumping occurs on mass scales of order 1 Mq. 
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1. INTRODUCTION 

Methanol is an abundant interstellar molecule, especially in regions of star formation where its 
high abundance is thought to be the result of thermal evaporation of dust grain mantles (Casein, 
Hasegawa, & Herbst 1993). A slightly asymmetric rotor with hindered internal rotation and sig- 
nificant a- and b-axis dipole moments, methanol has a complex spectrum which is sensitive to 
interstellar conditions. Observable lines cover a wide range of energies and line strengths, allowing 
for detailed analysis of excitation conditions and opacities. 

Two classes of interstellar methanol masers are known. Class I methanol masers, produced by 
collisional pumping, are seen in regions of massive star formation, but are generally well separated 
from the centers of activity as traced by embedded infrared sources and ultracompact H n (UC H n) 
regions. Class II methanol masers are thought to be produced by radiative excitation at infrared 
wavelengths and are spatially well correlated with young stellar objects, UCHii regions, and OH 
masers. Even when not masing, methanol may be expected to exhibit distinctly non-LTE excitation, 
except under the highest density conditions where collisions could thermalize the level populations. 

The southeastern portion of W3 contains the hydroxyl maser source W3(OH), which is associ- 
ated with an UC H n region around a young 07 star, and the nearby water maser source W3(H20), 
which contains a young stellar object known as the TW object (Turner & Welch 1984). Numer- 
ous studies have been made of molecular material associated with these two compact objects and 
throughout the region surrounding them (Wilson, Johnston, & Mauersberger 1991; Wilson et al. 
1993; Wink et al. 1994; Helmich k van Dishoeck 1997; Wyrowski et al. 1999; Nilsson et al. 2000). 
W3 is thought to be about 2.2 kpc from the Sun (Humphreys 1978). At such a distance, resolution 
of order 1 arcsec is necessary to study the nature of protostellar cores on scales of order 1000 AU. 

Sutton et al. (2001) presented results on new class II methanol masers in W3(OH) and discussed 
excitation conditions needed to explain observed fluxes in those and several other methanol maser- 
candidate lines. Sobolev, Sutton, & Watson (2003) extended that analysis to include absorption and 
emission components in a larger set of methanol lines in W3(OH) and discussed the implications in 
terms of excitation and source structure. Here we wish to complete our analysis of methanol in the 
vicinity of W3(OH) and W3(H2 0) by presenting our complete set of millimeter- wave observations 
and by discussing the nature of the methanol associated with the W3(H20) source and methanol 
in other nearby regions. 
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2. OBSERVATIONS 

Observations were carried out with the BIMA 1 interferometer between 1997 October and 2001 
June, following the general procedures discussed in Sutton et al. (2001). The set of observed 
12 CH30H transitions is presented in Table 1. In addition, observations were made of the central 
four lines of the J = 2-1 band of 13 CH30H. At 3 millimeters the data were primarily from the 
BIMA B and C arrays, with a small amount of A array data for two of the lines. At 1 millimeter the 
data were from the BIMA C and D arrays. In all cases spatial resolution was of order 2 arcsec, as 
shown in Table 1. Spectral resolution was generally of order 0.2 km s _1 , although some observations 
were made with resolution as wide as 1.2 km s -1 . The continuum was subtracted from the uv-data 
before maps were made. The phase center of the maps is a(J2000) = 2 h 27 m 03?87, (5(J2000) = 
61°52 / 24'. / 6. The flux density scale is thought to be accurate to within about ± 10% for the 3 mm 
band and about ± 25% for the 1 mm band. 

In order to check for spatially extended emission, the 96.7 GHz methanol J = 2-1 band also 
was mapped using the 20 meter Onsala telescope during 1998 January. The BIMA J = 2-1 maps 
and the Onsala maps are in general agreement with regards to the extent and shape of the emitting 
region. Furthermore, the BIMA interferometer maps recover about 75% of the peak flux density 
of these lines in the Onsala data. These lines are among the lowest energy lines in our data set, 
and therefore more likely to be excited in diffuse, spatially extended material. Since we appear to 
have recovered most of the flux in these low energy lines, we believe that we likely have done so 
also for the remaining lines in our data set. Since the BIMA maps have considerably higher spatial 
resolution and positional accuracy and since we only have Onsala data for four of our 38 lines, in 
the discussion which follows we treat the BIMA maps by themselves, without further reference to 
the Onsala data. 



3. RESULTS 

3.1. Low energy lines 

In Figure 1 we show maps from those methanol lines we have measured which have upper state 
energies less than 100 K and for which we have u-v data which are substantially complete down to 
spatial frequencies of 3 kilolambda or less. In other words, these maps are based on data sufficient 
to show any large scale structure which might be present. Diffuse emission is present from spatially 
extended material to the southwest of W3(OH). This southwest region will be discussed in detail 
in section 4.1, below. 



1 The Berkeley-Illinois-Maryland Association operates the BIMA array under support from the National Science 
Foundation. 
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3.2. Other a-type lines 

In addition to the four J = 2-1 lines from the ground torsional state, shown in Figure la,c,d, 
we also observed all six lines from the J = 2-1 band in the first torsionally excited state (vt=l). 
The torsionally excited lines are shown in Figure 2. These lines have upper state energies of order 
300-400 K and highlight the high excitation material in the vicinity of the two main sources of 
luminosity, W3(OH) and W3(H20). Taken together, these ten a-type lines from the V(=0,1 states 
form a standard against which other methanol lines can be compared, since the a-type lines are 
less affected by non-LTE conditions than the 28 b-type lines discussed below. 

We also observed the 2 -l A+, 2_i-l_i E, 2 -l E, and 2i-li E lines of 13 CH 3 OH. All four 
lines are weak, but visible towards both W3(OH) and W3(H 2 0). Flux densities are consistent with 
12 C/ 13 C isotopic ratios of order 60 and moderate opacities in the 12 CH30H lines, as discussed by 
A. M. Sobolev et al. (in preparation) and in section 4.2, below. 

3.3. Other b-type lines 

In Figure 3 we present maps of 15 b-type methanol lines, in addition to the 3 b-type lines 
already presented in Figure lb,e,f. Many of these lines are expected to be strong under class II 
(radiative pumping) conditions. These include the maser lines discussed in Sutton et al. (2001) 
(Fig. 3a,f,g) and the maser-candidate lines (Fig. lb; Fig. 3c,e,h,i,l) discussed in the same paper. 
Note the tendency for the emission to be stronger in these lines towards W3(OH) than towards 
W3(H 2 0), with the exceptions of the 6_ 2 -7_i E line at 85.5 GHz and the 8 3 -9 2 E line at 94.5 
GHz. Both of these lines are predicted to mase under some class II conditions by Sobolev, Cragg, 
& Godfrey (1997a), and the 85.5 GHz line has been seen masing in G345.01+1.79 (Sobolev et al. 
1998). However both lines were predicted to be in absorption under the specific conditions of high 
kinetic temperature derived by Sutton et al. (2001) for W3(OH). 

Figure 3 also contains data on three new transitions expected to be strong under radiative 
pumping conditions. The 5i~4 2 E line (Fig. 3b) has recently been seen as a class II maser in 
G345.01+1.79 by Sobolev et al. (2002), and indeed it is quite strong here in the vicinity of W3(OH), 
although there is no evidence of a narrow maser spike at —43.1 km s _1 , as was characteristic of 
the 86 and 107 GHz masers discussed by Sutton et al. (2001). However, the W3(OH) emission 
extending from about —48 to —43 km s" 1 could be a superposition of a number of weak maser 
features, similar to the case of several maser-candidate lines discussed in Sutton et al. (2001). 

The other new lines expected to be strong under radiative pumping are torsionally excited 
lines, lo~2i E vj=l and 6i-7 2 A - v<=l. These are particularly important in understanding the 
class II maser excitation mechanism, since that mechanism involves radiative excitation through the 
first two torsionally excited bands of methanol, v t =l and 2. Voronkov, Austin, & Sobolev (2002) 
recently observed the 2q-3i E vt=l line at 44.9 GHz in several galactic sources and concluded that 



- 5 - 



it is likely to be weakly masing in W3(0H). They predicted that the closely related lo~2i E v t =l 
93.1 GHz line would also be inverted in W3(OH), and our observations (Fig. 3m) demonstrate 
strong emission in that direction. Calculations show that 6i— 72 A - V(=l also is highly excited 
by the class II mechanism and can be a weak maser (A. M. Sobolev et al., in preparation). We 
do not have enough evidence here to confirm masing in either of these lines. But the calculations 
mentioned above indicate super-thermal excitation in both lines in the vicinity of W3(OH). 

Other b-type lines are expected to be favored under class I (collisional pumping) conditions. 
The 5_i-4o E and 8o~7i A + lines will be discussed extensively in §4.1. Here in Figure le,f we draw 
attention to the fact that emission from these lines is somewhat stronger from W3(H20) than from 
W3(OH). Data are presented on three other lines from the ground torsional state: 8o~7i E, 11 i 
10_2 E, and 13i-122 A + (Figures 3d, j, and k, respectively). It is unclear how strongly 8o— 7i E and 
13i-122 A + might be favored by class I conditions and disfavored by class II conditions, although 
the data show that 13i-122 A + emission is indeed stronger towards W3(H20) than W3(OH). The 
104.3 GHz ll_i~10_2 E line belongs to the same series as the 9.9 GHz class I methanol maser 
9_i-8_2 E. To the best of our knowledge, the 9.9 GHz maser has been seen, weakly, in only one 
source, W33-Met (Slysh, Kalenskii, & Val'tts 1993). Voronkov and Alakoz have recently seen maser 
emission in the 104.3 GHz line, again only in W33-Met (M. Voronkov 2003, private communication). 
Calculations reported by Cragg et al. (2001) suggest that the relative excitation of the K = — 1 and 
K = — 2 E ladders is very model dependent. Here we simply note that emission from ll_i-10_2 E is 
stronger towards W3(OH) than W3(H20). Finally in Figure 3n we present data from a torsionally 
excited line, 6i~5o E vj=l, which should be more favored by class I conditions than class II, and 
indeed it is stronger towards W3(H20) than W3(OH). 

We have data on ten other b-type lines, listed in Table 1 but not presented in Figures 1 or 
3. Although these additional lines will be used in the analysis which follows, the maps either 
have relatively poor sensitivity or the lines are not strongly preferred by either class I or class II 
conditions. For example, four of these lines are relatively high J members of the J_2-Ji E species 
Q-branch. Most of the remaining lines in this category are of similarly high J or have |K| > 3. 
They are of minor importance in the analysis. 

4. DISCUSSION 

4.1. Spatially extended material to the southwest 

Some emission to the southwest of W3(OH) is visible in all of the lines shown in Figure 1, 
although it is most extensive in the lowest energy lines plotted in Figure la. This is a blend of the 
2o~lo A + (E u = 7 K) and 2_i-l_i E (E u = 12.5 K) lines. The relative weakness of the emission in 
this region for the other lines of the J = 2-1 band, 2 -lo E at E u = 20 K (Fig. lc) and 2i-li E at E u 
= 28 K (Fig. Id), indicates this is optically thin emission from low excitation material. An average 
of the emission in the lower right quadrant of the plots for the J = 2-1 lines (excluding emission 
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within 4 arcsec of W3(0H) itself) implies a rotational temperature of 11 K. This is consistent with a 
kinetic temperature of order 50 K and a density of order 10 5 -10 6 cm~ 3 (more detailed calculations 
are given below). This methanol emission shows good spatial and kinematical correspondence with 
material to the southwest of W3(OH) seen in NH3 (1,1) and (2,2) emission by Wilson et al. (1993), 
who find the gas in that region to be characterized by a hydrogen density of 6 x 10 5 cm -3 and a 
rotational temperature of 30 K. It appears that in methanol we are seeing material similar to that 
seen in metastable ammonia. The presence of high density material in this region is also suggested 
by the SO and CS data of Nilsson, Bergman, & Hjalmarson (2000). 

Within the low excitation methanol region are two clumps which show up prominently in the 
5_i-4o E and 8o-7i A + lines. These are the clumps labelled swA and swB in Fig. le,f, at offsets 
of (Aa = -8, AS = -14) and (Aa = -13, AS = -7) arcsec, respectively. Since W3(OH) itself is a 
class II methanol maser source, the presence of strong compact emission in these lines, known for 
their strong class I maser emission, is unexpected. It seems that in these two regions at least, the 
excitation conditions are more like those found in class I sources. The class I mechanism is driven 
by collisional pumping. Therefore, in the regions of clumps swA and swB, collisions evidently are 
producing super-thermal excitation in these, and potentially other, class I methanol maser lines, 
whereas radiative pumping is less important. These positions are also near local maxima in the 
2o-lo A + and 2_i-l_i E map (Fig. la). The class I mechanism can also modestly enhance the 
intensity of these two particular J = 2-1 lines, although it seems likely that these positions are also 
local maxima in methanol column density. 

To show that either of these regions is actually weakly masing via the class I mechanism 
requires further evidence such as sub-thermal linewidths or brightness temperatures well in excess 
of conceivable kinetic temperatures. Data on linewidths and flux densities from clumps swA and 
swB are contained in Table 2. The linewidths are narrow (2-3 km s~ x FWHM for the 5_i~4o E 
and 80— 7i A + lines), but not sufficiently so to conclude that these are masers (linewidths of class 
II masers in this source are typically 0.3 km s _1 , and thermal widths at 50 K are of similar order). 

Emission is strong and compact for clump swA in the 80— 7 1 A + line, so this would seem to be 
the best case for making an argument in favor of maser emission based on brightness. The spatial 
resolution of the map show in Figure If is 2.5 x 2.2 arcsec. Clump swA is partially resolved in this 
map. Based on a deconvolved source size of 3.0 x 1.9 arcsec and a peak flux density of 2.0 Jy, we 
derive a brightness temperature of 48 K. If structure is present at spatial scales smaller than the 
current resolution, this represents a lower limit on the brightness. This is not sufficient to prove 
class I maser action, but since this temperature is greater than that known for other material to 
the southwest of W3(OH) and there is no known source of luminosity near this position, it seems 
likely that some degree of maser amplification is taking place. The alternative, that this is simply 
a high density clump in otherwise rather warm material, seems unlikely due to the weakness of the 
2o-lo E and 2\-li E emission in this region. 

The remaining line to be discussed in this context is 0q-1-i E, presented as Figure lb. In 
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this transition we see relatively little emission at the locations of clumps swA and swB as well as 
throughout the region to the southwest of W3(OH), even though the upper state energy of this 
transition is only 13.1 K. The reason for this appears to be that this line is favored by class II 
excitation conditions and is a known class II maser in G345.01+1.79 (Val'tts et al. 1999). Lines 
which mase or are super-thermally excited under class II conditions are expected to be sub-thermally 
excited under class I conditions such as prevail here to the southwest of W3(OH) and specifically 
near clumps swA and swB. 

It might seem difficult to be precise about the conditions in these regions due to the small 
number of lines detected towards clumps swA and swB, the likely non-thermal nature of the exci- 
tation, and the possibility of significant opacity in the lowest energy lines. However, non-LTE fits to 
the data for these seven lines do provide significant constraints on source parameters, as shown in 
Table 3. These calculations were performed according to the methods described by Sobolev, Cragg, 
& Godfrey (1997a), Sutton et al. (2001), and Salii, Sobolev, & Kalinina (2002), but without an 
incident infrared radiation field. Assuming the warm dust is localized to the immediate vicinities 
of W3(OH) and W3(H20), the infrared field should be down by a factor of at least 500 at the 
positions of swA and swB, relative to that near W3(OH). The physical parameters of kinetic tem- 
perature, hydrogen density, methanol column density, and source size were varied to obtain optimal 
agreement between the model and the observed fluxes. Parameter uncertainties were obtained in 
the usual fashion (Press et al. 1992) from the curvature of the x 2 hypersurface. Each represents a 
68% confidence limit after marginalizing over the remaining parameters. 

Comparing our methanol column densities with an estimated hydrogen column density of 
3 x 10 23 cm~ 2 for this region (Wilson et al. 1993), we derive a methanol fractional abundance of 
order 10~ 8 . Independent estimates of hydrogen column densities for the individual sources can be 
obtained by multiplying the hydrogen number densities from Table 3 by the thicknesses along the 
lines of sight (assumed to be equal to the derived source diameters of 4400 AU for swA and 5000 
AU for swB). These yield methanol fractional abundances of 3.7 x 10~ 7 for swA and 3.1 x 10 -8 for 
swB. These methanol fractional abundances are higher than the values of order 10~ 9 usually seen 
in dark molecular clouds (Kalenskii & Sobolev 1994). But they appear to be lower than the values 
of order 10 -6 found for other star forming regions (Menten et al. 1986) and generally attributed 
to thermal evaporation of dust grain mantles. Thermal evaporation of methanol from polar ice 
mantles is thought to occur at temperatures of 85 - 90 K (Hasegawa &: Herbst 1993; van Dishoeck 
et al. 1993; Tielens & Whittet 1997; van Dishoeck & Blake 1998). Therefore the higher methanol 
fractional abundance in swA is consistent with its higher temperature, as shown in Table 3. At 
the minimum projected distances of these clumps from W3(OH), neither gas nor dust grains are 
likely to reach temperatures of this magnitude due solely to the luminosity of W3(OH) (Doty & 
Neufeld 1997). Therefore it is possible that swA contains an internal heating source. Since the 
temperature of swB is insufficient to evaporate methanol, some non-thermal desorption mechanism 
may be necessary to release the observed amount of methanol. 

Assuming uniform densities, masses for the gas in clumps swA and swB also may be obtained, 
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as shown in Table 3. These masses are of order 0.05 M Q and 1.1 M , respectively, with large 
uncertainties which are probably dominated by uncertainties in the sizes and filling factors. These 
numbers may be compared with an estimate for the total mass in a larger region to the southwest 
of W3(OH) of 80 M Q (Wilson et al. 1993). The observed linewidths imply virial masses of order 
25 M each for swA and swB. It is unlikely that either clump is gravitationally bound, unless it 
contains a compact core of significant mass which remains unresolved in the current observations. A 
collapsing high mass core or a high mass protostar at a very early evolutionary stage, in the center 
of swA, could provide sufficient mass to bind the surrounding gas cloud as well as the luminosity 
needed to produce the observed gas temperature. 

For source swA, our best model fit shows that all of the lines listed in Table 2 are inverted, 
except for Oo~l-i E. However, only 5_i-4o E and 8o~ 7i A + show significant amplification, with 
r = —0.8 and r = —0.9, respectively. For source swB the energy levels are more nearly thermalized, 
as expected for the higher H2 density. Although three of the lines are still inverted, only 8o-7i 
A + shows significant amplification, with r = —0.6. Our model predicts that at these locations 
the closely related 36 GHz (4_i-3o E) and 44 GHz (7o~6i A + ) methanol maser lines will also be 
inverted and show modest amplification (r ~ —1)- Higher J members of these series exhibit weaker 
inversions. Observed line widths of 2.1 km s _1 (Haschick & Baan 1989, 36 GHz), 0.8 km s _1 
(Haschick, Menten, & Baan 1990, 44 GHz), 2.1 km s" 1 (for swA at 84.5 GHz; this work), and 2.5 
km s _1 (for swA at 95.2 GHz; this work) are narrower than typical widths for this source (Kalenskii 
et al. 1997), consistent with modest amplification in these lines. 

Our modelling of the southwest region gives a possible explanation for some of the complex 
behavior of the 25 GHz J2-J1 E lines near W3(OH). The observations of Menten et al. (1986) showed 
a mixture of absorption and emission line components in this band when viewed with a large beam 
towards W3(OH). The strongest emission was seen for J=3, although emission was evident up to 
J=9. Our model of swB predicts a weak inversion in the J=3 line, super-thermal excitation in 
J=4, and weak emission for other nearby values of J. In contrast, Sutton et al. (2001) predicted 
absorption from the material directly in front of the W3(OH) UCH11 region. Since the Menten et 
al. (1986) beam encompassed both W3(OH) and swB, their observations showed a mixture of these 
characteristics, further complicated by emission and absorption from other methanol-rich regions 
contained within their beam (cf. Fig. la). 

4.2. W3(H 2 0) region 

Most of the lines studied here exhibit strong emission from the region around the water maser 
source W3(H20). Careful examination of our maps with the best signal-to-noise ratio and highest 
spatial resolution reveal that the bulk of the methanol emission comes from a region approximately 
2.0 x 1.2 arcsec in size, elongated along a position angle of 100° ± 10° and centered about 5.3 arcsec 
east of our nominal map center (Aa = 5'/3 ± O'/l, A5 = O'/l ± O'/l). In other words, the emission 
comes from a region encompassing the TW object (Turner k, Welch 1984) and continuum peaks B 
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& C of Wyrowski et al. (1999). 

This is illustrated in Figure 4, which presents a 0.6 arcsec resolution map of the 1.3 mm 
continuum from this region. Sources A (the TW object), B, and C of Wyrowski et al. (1999) are 
labelled. Our continuum map has slightly lower resolution than Wyrowski et al., and we do not 
resolve their sources B and C from each other. A solid line is overlaid on the map to represent the 
elliptical region of methanol emission described above. The extent of the methanol emission roughly 
corresponds to the 0.2 Jy beam -1 contour of the continuum map. The similarity in position, size, 
shape, and orientation of the methanol and the continuum emission strongly suggests they are 
coming from the same material. We note that the agreement in position can be made even better 
by shifting the 1.3 mm continuum map south 0.1 arcsec and west 0.1 arcsec. Such a shift is within 
the estimated positioning error of our map and also serves to shift the center of the W3(OH) 1 
mm continuum onto our previously reported 3 mm position, as shown (Sutton et al. 2001). It also 
brings our two W3(H20) peaks into better agreement with peak A and an average of peaks B and 
C of Wyrowski et al. (1999). 

In Table 4 we present the observed integrated line fluxes for these lines, and in Figure 5a we 
plot these fluxes versus upper state energies in the usual "rotation diagram" format. The diagram 
assumes a source solid angle of $7 = 4.4 x 10~ n ster, a value which gives the optimal fit in the 
procedure described below. In the limit of optically thin emission from a region in LTE, the points 
in a rotation diagram should fall on a straight line with the slope indicating the kinetic temperature 
and the intercept containing information about the total column density, according to the relations 

N u 3c 2 1 / S^dv 

~ i6vr-v 3 n fi 2 s 

an d 

N„_ = N T e „ Eu/kTrot 
gu Q(T rot ) 

In these equations the total column density, Nt, and the upper state column density divided by 
the level degeneracy, N u /g u , are averages over the source solid angle, Q. The rotational partition 
function is Q(T rot ), and the vibrational partition function is neglected. The observed flux integrated 
over velocity is f S^dv. The energy of the upper state is E u , and the square of the dipole moment 
matrix element for the transition is fi 2 S. Values for both are given in Table 1. The other symbols 
take on their usual meanings. The data show large deviations from straight line behavior, indicating 
the breakdown of one or more of the assumptions implicit in the rotation diagram technique. Several 
of the lines shown are in fact optically thick. In a rotation diagram, optically thick lines appear 
artificially low since the observed intensity is less than that expected from a linear extrapolation 
from the optically thin regime. The optical depth of a line in LTE is given by 

T - — // 2 S ^ p -E u /kT rot / hv/kT Iot _ 1 x J_ 

T " 3h M b g 1 ' Av ' 

where Av is the linewidth and No/go is the column density of the ground state. The optical 
depths of these lines range from 0.1 to 4.6, with approximately half of the lines showing significant 
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saturation. The effects of saturation may be corrected using the factor C T (Goldsmith k, Langer 
1999; Gibb et al. 2000), where 



In Figure 5b we present a "corrected" rotation diagram, in which the deviations from straight line 
behavior are much reduced. It is apparent that the largest deviations in the original diagram were 
the result of saturation. The best fit in Figure 5b has an intercept of No/go = 1-3 x 10 15 cm~ 2 , a 
slope corresponding to a rotational temperature of 149 K, and a solid angle of 4.4 x 10~ n ster. This 
gives a value of x 2 = 128 for 35 degrees of freedom (38 data points minus three free parameters). 
Some significant deviations from a good fit still remain, so the next step is to determine whether 
they represent real departures from LTE. 

A non-LTE environment in which collisional pumping dominates over radiative excitation 
produces conditions similar to those in class I methanol masers. Under these conditions, certain 
lines are predicted to have super-thermal excitation, particularly those with upper states in the 
K = — 1 E and K = A ladders. Four out of the five such lines in this data set lie above the 
average trend of the data in Figure 5b. Similarly, b-type lines with lower states in these ladders are 
predicted to show sub-thermal excitation. Again four out of five such lines lie below the average 
trend of the data. Therefore, the conditions around W3(H20) appear similar to those in class I 
maser regions, even though W3(H20) is not a methanol maser source. 

We performed detailed non-LTE calculations for the data in Table 4, using a simplified model 
of the infrared radiation field. Such a treatment is justified since at these densities methanol is near 
the collisionally dominated limit. Ostrovskii & Sobolev (2002) have shown that when the color 
temperature of the radiation is in the range 125-250 K, the methanol radiative excitation occurs 
primarily through a large number of strong methanol absorption lines in the 15-40 //m wavelength 
range. The most important parameter of the infrared radiation field therefore is its intensity in 
this band. We assume that the pumping radiation is produced by warm dust intermixed with the 
gas and at the same physical temperature as the gas. At these wavelengths the dust is optically 
thick in this source, minimizing the importance of infrared sources external to the cloud (although 
the TW source and Wyrowski's B/C remain the ultimate sources of energy in this region). 

Methanol radiative transition rates were taken from Mekhtiev, Godfrey, & Hougen (1999). 
Collisional transition rates were based on the model of Peng & Whiteoak (1993). For simplicity, 
a A~ 2 dust opacity law was chosen. Details of the excitation do depend on the optical properties 
of the dust grains (Ostrovskii & Sobolev 2002), but the differences are small at high gas densities. 
Radiative transfer for the millimeter-wavelength methanol lines was treated using the large velocity 
gradient (LVG) approximation. Further details regarding the calculational methods may be found 
in K6ppen & Kegel (1980), Sobolev (1982), Sobolev, Cragg, & Godfrey (1997a), and Sobolev, 
Cragg, & Godfrey (1997b). 

The gas cloud was assumed to be uniform in H2 density, gas and dust temperature, dust- 
to-gas ratio, and methanol fractional abundance. The data were not adequate to justify a more 
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complicated, spatially inhomogeneous model of the cloud. The source solid angle was fixed at 
4.4 x 10" 11 ster, as given above. The other four parameters were allowed to vary. An optimal 
fit for W3(H2 0) was obtained with a kinetic temperature of 140 ± 10 K, hydrogen density of 
0.9±l$ x 10 7 cm" 3 , and N A + N E = 2.4+°; 3 x 10 18 cm" 2 , as shown in Table 3. We can multiply 
the hydrogen number density by the estimated source thickness of 3400 AU to derive a hydrogen 
column density. Combining this with the measured methanol column density gives a methanol 
fractional abundance of order 5 x 10~ 6 for W3(H20). In comparison with the values in Table 3, 
van der Tak et al. (2000a) give a kinetic temperature of greater than 200 K for W3(H20) from 
statistical equilibrium calculations for H2CO, similar to those given earlier by Helmich et al. (1994) 
and Helmich & van Dishoeck (1997). However, the beam size for those observations was large 
enough to encompass W3(OH) as well as W3(H20). The methanol column density of van der Tak, 
van Dishoeck, & Casern (2000b) is smaller than that given in Table 3, an effect we also attribute 
to their larger beam size. The density profile of van der Tak et al. (2000a) corresponds to a density 
of about 3 x 10 7 cm" 3 at a characteristic radius of 2000 AU, consistent with our result. 

The observations also provided good constraints on the amount of dust needed in the model. 
This gives us an independent method to calculate the hydrogen column density and derive a 
methanol fractional abundance. Previously we adopted a A" 2 dust opacity law. Now, only for 
the purpose of this calculation, we further adopt a gas-to-dust mass ratio of 100 and a dust mass 
absorption coefficient of 1.1 cm 2 g" 1 at 1.3 mm (Henning, Michel, & Stognienko 1995). This yields 
a methanol fractional abundance of order 1.7 x 10" 6 for W3(H20). Given the uncertainties inherent 
in these various assumptions, the two estimates of the fractional abundance agree well and support 
the notion of significant grain evaporation in this source. 

Figure 5c shows the residuals from the optimal non-LTE fit. For this fit x 2 = 86 for 34 degrees 
of freedom (38 data points minus 4 free parameters). The non-LTE model is a significantly better 
fit to the data than the LTE model. There are still some significant residuals, possibly indicating 
that the observational uncertainties were underestimated. But it should be kept in mind that 
there are additional uncertainties, those associated with the model calculations. Such uncertainties 
are difficult to quantify, but they appear to be of the same order as those associated with the 
observations. The nature of these uncertainties is discussed in Appendix B of Sutton et al. (2001). 
The dominant uncertainty may be spatial variation in the source parameters. 

The results presented above are based on averages over the entire W3(H20) source. Although 
our spatial resolution is generally insufficient to resolve the TW object from sources B and C of 
Wyrowski et al. (1999), we do have some information on spatial sub-structure from the shape of 
the spectra. For the 29 v^=0 lines listed in Table 1, the average methanol emission velocity towards 
W3(H20) is —50.0 km s -1 , with a standard deviation of 0.37 km s" 1 (giving a standard error of 
the mean of less than 0.1 km s" 1 ). Several lines with good signal-to-noise ratios show distinct 
shifts from this velocity, as shown in Figure 6. The 6i-72 A™ v t =l line is well fit by a single 
gaussian with a velocity of — 51.0±0.1 km s -1 . This blueshifted emission is more characteristic 
of the TW source in the eastern part of W3(H20). Note that in Figure 3o the emission peak for 
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this line is spatially coincident with the TW object. This is a high excitation line (E u = 374 K), 
presumably more easily excited close to the main source of luminosity in the region, namely the 
TW object. The 8o~7i A + and 5_i-4o E lines have mean velocities of — 49.3±0.1 and — 49.5=L0.1 
km s" 1 , respectively. Redshifted emission is characteristic of the western portion of W3(H20), near 
Wyrowski's source C. And in the maps of these lines, the emission peaks distinctly to the west of 
TW (Figures le,f - although this shift to the west is too small to be seen on the scale at which these 
maps are presented). So the data suggest that the line-of-sight velocity difference between the TW 
object and source C may be of order 1.6 km s _1 . The 8o~7i A + and 5_i— 4o E lines are the same 
two lines discussed in section 4.1, known elsewhere for their class I methanol maser emission. Since 
such lines are known to be collisionally pumped, the luminosity source proposed by Wyrowski et al. 
(1999) may either be sufficiently weak or sufficiently well shielded from the methanol-rich material 
for collisional pumping to dominate. We should note that in our non-LTE model for the W3(H20) 
region, these two lines are optically thick, with opacities of ~3.7 and ~1.6, respectively. So we may 
be seeing a preponderance of foreground material in these lines. 

Given a density of 0.9 x 10 7 cm -3 and assuming that W3(H20) is as deep along the line 
of sight as its average extent in the plane of the sky (3400 AU), this corresponds to 1.6 M Q of 
warm gas, excluding that contained in proto-stellar sized objects and condensed cores. This may 
be compared with a previous estimate of the gas mass of "at least a few M " (Turner & Welch 
1984). The estimate of 70 M by Wilson, Johnston, & Mauersberger (1991) refers to a considerably 
larger region (0.13 pc) probed by C 18 0. Wink et al. (1994) estimate 60 M Q in a 0.05 pc diameter 
region from C 18 observations; however, they also estimate 30 M on a 0.01 pc scale from dust 
observations. Finally, Wilson et al. (1993) give a dynamical mass of 18 M , which therefore includes 
that of any protostars. 

Non-thermal continuum emission at centimeter and decimeter wavelengths has been seen both 
east and west of the TW object (Reid et al. 1995; Wilner, Reid, Sz Menten 1999; Argon, Reid, & 
Menten 2003). This has been interpreted as synchrotron radiation from precessing jets (Wilner, 
Reid, & Menten 1999), in part due to the symmetry and curvature of the emission and its location 
near the apparent center of expansion of the H2O masers. However, alternative explanations are 
possible. Torrefies et al. (2002) point out that in some cases of star formation, water masers are 
found in protoplanetary accretion disks, oriented perpendicular to the radio jet directions (although 
Torrefies et al. (2002) do not consider the TW object to be such a system). But recently Shchekinov 
&; Sobolev (2003) have proposed a model in which the H2O masers in W3(H20) are formed on the 
surface of an east- west circumstellar disk or ring. 

This question of orientation has some bearing on the nature of the object we have been calling 
Wyrowski's source C. Its location is along the east-west axis of of the synchrotron jet proposed 
by Reid et al. (1995). In fact, it is roughly coincident with the cluster of water masers located 1 
arcsec west of the TW object. Since jets produce shocks where they run into ambient material, and 
shocks, in turn, can power masers, one possibility is that Wyrowski's source C is a deeply embedded 
Herbig-Haro object. However, Wilner, Reid, & Menten (1999) and Wyrowski et al. (1999) agree 
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that regardless of jet geometry, the most likely interpretation is that Wyrowski's source C is a 
deeply embedded protostellar object. We concur. 

Using this interpretation and assuming that the molecular gas seen in methanol is roughly 
equally divided between the TW object and Wyrowski's source C, we get the following picture. The 
TW object itself is a protostar of about 15 M Q (Mauersberger, Wilson, & Henkel 1988). Wyrowski's 
sources B/C may contain additional protostellar cores of similar or smaller mass (Wyrowski et al. 
1999). Within 0.01 pc of each of these sources there is about one solar mass of warm gas. And on 
larger scales there are an additional few tens of solar masses of cooler, less dense gas. 

The dynamical stability of such a system is uncertain. The putative sum of the protostellar 
masses (20-30 M ), the projected separation («2500 AU), and the estimated line of sight veloc- 
ity difference (1.6 km s _1 ) are consistent with a reasonable binary orbit of the TW object with 
Wyrowski's B/C. For example, a low eccentricity, high inclination, intermediate phase (an angle of 
order 45° with respect to the line of nodes) orbit is consistent with the above values. This assumes 
that the 1.6 km s _1 velocity difference discussed above in connection with Figure 6 represents the 
true difference in line of sight protostellar velocities. This would be the case if each protostellar 
object were accompanied by its own circumstellar material, but not necessarily the case if they were 
orbiting within a common, circumbinary envelope. A numerical simulation by Turner et al. (1995) 
produces a binary system with circumstellar material in the form of individual protostellar disks, 
a possible case of the system described above. Lubow &; Ogilvie (2000) indicate that individual 
disks of sufficient size, like those of Turner et al. (1995), can be stable in binary systems. The full 
width to zero power linewidths seen here are 10-12 km s^ 1 . Such velocities could be produced by 
slow infalls or outflows or by orbital motion in individual circumstellar envelopes. A significant 
amount of circumbinary material is also possible, although if it is in purely orbital motion it could 
not account for the full widths of these lines. 

In a single dish study of high mass star forming regions, Minier & Booth (2002) report emission 
from the vicinity of W3(OH/H20) in a number of methanol lines. They report emission peaks at 
-51 km s" 1 in the 6_ 2 -7_i E (85.5 GHz) and 7 2 -8i A+ (111.2 GHz) lines. Our work shows that 
gas at this velocity is clearly associated with the W3(H20) source. Minier & Booth (2002) also 
point out that W3(OH/H2 0) is unusual among massive star forming regions in not displaying clear 
maser emission in the 8o~7i A + (95.1 GHz) line. The emission in this line from sources swA and 
swB was discussed above in section 4.1 and that from W3(OH) by Sobolev, Sutton, & Watson 
(2003). For the W3(H20) source our model suggests an excitation temperature for this line of 
about 135 K and an optical depth of about 3.7, consistent with the lack of maser emission. 

In our optimal non-LTE model of W3(H20) none of the lines studied here are inverted; exci- 
tation temperatures range from about 90-530 K. Optical depths of these lines range from about 
0.07 to 3.9. 

In our optimal non-LTE model, the lines in the J = 2-1 band of 12 CH30H have opacities 
ranging from 1 to 1.5. This means that, even apart from non-LTE effects, the flux ratios for 
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the J = 2-1 band of CH3OH relative to 3 CH30H would be expected to be smaller than the 
actual abundance ratio by a factor of 1.9. In addition, since the 13 CH30H lines are optically 
thin, 13 CH30H will be more efficiently cooled, resulting in higher fractional populations in the 
lower rotational levels. A detailed calculation of this effect would require a larger set of 13 CH30H 
observations and a separate non-LTE model for 13 CH 3 OH. However, the cooling is likely to be 
substantial. The [ 12 C]/[ 13 C] abundance ratio is thought to be 66 ±4 (Langer & Penzias 1990). Due 
to the above effects, the integrated flux ratio for the J = 2-1 band should be considerably smaller, 
possibly consistent with the observed ratio of 12 ± 3, ignoring fractionation. 

5. Conclusions 

To the southwest of W3(OH) there is an extensive region of gas, dense enough to be seen in 
a number of methanol lines. Emission is most readily seen in lines of low upper state energy (E u 
< 15 K) and in lines known elsewhere for their class I methanol maser emission. In this case, 
strong class I masers do not seem to be present, but the physical conditions are similar to those 
seen in class I maser regions. Excitation is collisionally driven, with densities of order 10 5 - 10 6 
cm -3 . Inverted populations and weak amplification appear to be present. We predict that the 
36 GHz and 44 GHz lines will also be inverted, although it is doubtful that their optical depths 
will be sufficient to produce strong masers. The region is clumpy on scales of about 1 to a few 
tens of solar masses. Kinetic temperatures range up to about 100 K. The warmest regions have 
methanol fractional abundances of order 10 -7 , suggesting that some grain mantle evaporation has 
taken place, but that the evaporation is incomplete. 

In contrast, towards W3(H20) the methanol emission is nearly thermalized and many lines 
are optically thick. None of the lines studied have inverted populations. Excitation is due to a 
combination of collisional and radiative processes. The size of the methanol emitting region is 
about 2600 x 4400 AU. This is resolved into at least two components, corresponding to the TW 
object and a blend of Wyrowski's B and C components. It is likely that both are massive young 
stellar objects. Their separation and line of sight velocity difference are consistent with a wide 
binary orbit with a period of order 50,000 years. Each source may be accompanied by circumstellar 
material containing of order one solar mass or less of gas. There may be as much as a few tens of 
Mq of circumbinary material, although its distribution is unclear. Methanol fractional abundance 
in the circumstellar material is of order 10~ 6 , indicative of a warm environment and extensive grain 
mantle evaporation. 
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were supported by INTAS grant 99-1667, RFBR grant 03-02-16433, and the Russian Ministry of 
Education (grant E02-11.0-43). The authors thank Yuri Shchekinov for helpful discussions. 
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Table 1. Observed 12 CH 3 OH Lines 



Transition 




Frequency 




ti 2 S 


Beam 






(MHz) 


(K) 


(D 2 ) 


(arcsec 2 ) 


7 2 -8i A - 




80993.26 


103. 


2.52 


2.6 x 2.5 


13-3-14-2 E 




84423.71 


274. 


4.30 


2.5 x 2.4 


5_i-4 E 




84521.21 


40. 


3.08 


2.5 x 2.4 


6_2-7-i E 




85568.07 


75. 


2.01 


3.7 x 2.7 


7 2 -6 3 A" 




86615.60 


103. 


1.36 


2.0 x 1.7 


7 2 -6 3 A+ 




86902.95 


103. 


1.36 


2.0 x 1.8 


8-4-9-3 E 




89505.78 


171. 


1.56 


3.6 x 2.6 


l -2i E v t = 


1 


93196.66 


303. 


1.34 


2.3 x 2.1 


83—92 E 




94541.81 


131. 


2.24 


2.8 x 2.2 


8 -7i A+ 




95169.52 


84. 


7.22 


2.5 x 2.2 


2i-li A+ v t = 


1 


96396.06 


332. 


1.21 


2.1 x 1.9 


2i-li E v t = 


1 


96492.16 


298. 


1.21 


2.1 x 1.9 


2 -lo E v t = 


1 


96493.55 


308. 


1.62 


2.1 x 1.9 


2_i-l_i E v t = 


1 


96501.70 


420. 


1.21 


2.1 x 1.9 


2 -lo A + v t = 


1 


96513.67 


431. 


1.62 


2.1 x 1.9 


2i-li A" v*= 


■1 


96588.59 


332. 


1.21 


2.3 x 2.0 


2_i-l_i E 




96739.36 


12.5 


1.21 


2.1 x 1.8 


2 -lo A+ 




96741.38 


7.0 


1.62 


2.1 x 1.8 


2o _ lo E 




96744.55 


20. 


1.62 


2.1 x 1.8 


2i-li E 




96755.51 


28. 


1.24 


2.1 x 1.8 


6i-5o E vt= 


1 


99730.96 


340. 


3.14 


2.1 x 1.8 


13 2 -12 3 E 




100638.87 


234. 


3.84 


3.2 x 2.4 


13-3-12-4 E 




104060.72 


274. 


3.27 


3.7 x 2.5 


II-1-IO-2 E 




104300.40 


159. 


3.41 


3.1 x 2.3 


13_ 2 -13i E 




104336.64 


237. 


1.20 


2.9 x 1.5 


IO4-H3 A" 




104354.86 


208. 


2.48 


3.1 x 2.3 


IO4-H3 A+ 




104410.49 


208. 


2.48 


3.1 x 2.3 


13i-12 2 A+ 




105063.76 


224. 


4.27 


2.3 x 1.9 


14_2-14i E 




105576.39 


270. 


1.79 


2.3 x 1.9 


3i-4 A+ 




107013.81 


28. 


3.01 


1.6 x 1.4 


15_2-15i E 




107159.92 


305. 


2.61 


1.8 x 1.6 


Oo-l-i E 




108893.93 


13.1 


0.98 


2.2 x 1.8 


14 5 -15 4 E 




109138.69 


380. 


3.41 


2.2 x 1.9 


I6-2-I61 E 




109153.21 


342. 


3.68 


2.2 x 1.8 


7 2 -8i A+ 




111289.60 


103. 


2.34 


1.9 x 1.6 


5i-4 2 E 




216945.56 


56. 


1.12 


2.5 x 2.3 


6i-7 2 A" v t = 


= 1 


217299.16 


374. 


4.66 


2.5 x 2.3 


80— 7 1 E 




220078.49 


97. 


3.45 


2.4 x 2.2 
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References. — Xu and Lovas 1997, except the 107103.81 
MHz line is from Tsunekawa et al. 1995. 
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Table 2 


. Data for clumps 


swA & swB 




Transition 


v lsr 


Av 


o 

°peak 


/ S„ dv 


Model 




(km s- 1 ) 


(km s- 1 ) 


(Jy) 


(Jy km s _1 ) 


(Jy km s _1 ) 




clump swA 


2 -l A+ 


-45.9 


3.3 


1.8 


6.0±0.6 


5.7 


2-i-l-i E 


-45.9 


3.6 


1.4 


5.1±0.5 


4.9 


Oo-l-i E 


-46.7 


4.6 


0.4 


1.6±0.3 


0.8 


2o-lo E 


-46.3 


3.4 


0.6 


2.0±0.2 


2.3 


2i-li E 


-46.3 


3.1 


0.3 


1.0±0.2 


0.4 


5_i-4 E 


-45.8 


2.1 


2.2 


5.3±0.5 


5.0 


8 -7i A+ 


-46.4 


2.5 


2.0 


5.2±0.5 


5.6 


clump swB 


2 -l A+ 


-46.4 


2.8 


3.1 


9.5±1.0 


8.1 


2-i-l-i E 


-46.6 


2.5 


2.5 


6.7±0.7 


6.5 


Oo-l-i E 


-46.6 


2.4 


1.4 


3.6±0.4 


3.2 


2o-lo E 


-46.4 


2.3 


1.8 


4.3±0.4 


5.1 


2i-li E 


-46.4 


2.3 


1.0 


2.4±0.3 


2.1 


5_i-4 E 


-46.4 


2.2 


4.0 


9.4±1.0 


4.5 


8 -7i A+ 


-46.2 


2.7 


2.2 


6.1±0.6 


6.7 
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Table 3. Physical parameters for clumps swA & swB and W3(H20) 



Source 


Tkin 

(K) 


n H 2 

(cm -3 ) 


N A +N E 
(cm -2 ) 


Mgas 

(M ) 


(N A +N E )/N H2 


swA 


103±34 


1.4l^ x 10 5 


3.4± 2 ; 3 x 10 15 


» 0.05 


3.71^ x 10- 7 a 


swB 


47± 6 


2.0lf J x 10 6 


4,.6t\i x 1015 


» 1.1 


3.lt L r x 10- 8 a 


W3(H 2 0) 


140±10 


0.9ig;| x 10 7 


2A±° i x 10 18 


» 1.6 


5. ±f x 10~ 6 a 

1.7l^ x io~ 6 b 



a N Ha is based on density and estimated path length. The uncertainty in the 
fractional abundance typically is dominated by the uncertainty in density. 

b Estimated fractional abundance based on infrared pumping. The stated un- 
certainty reflects only that associated with the amount of infrared pumping and 
ignores uncertainties in the gas-to-dust ratio, the dust absorption coefficient, and 
the dust opacity law. 

Note. — Estimated 68% confidence limits. 
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Table 4. Integrated Fluxes from W3(H 2 0) 



Transition 




Frequency 


En 


/ S„ dv (obs.) a 


Model 






(MHz) 


(K) 


(Jy km s _1 ) 


(Jy km s _1 ) 


7 2 -8i A" 




80993.26 


103. 


3.5± 0.4 


3.7 


13_3-14_2 E 




84423.71 


274. 


2.2± 0.2 


2.5 


5_i-4 E 




84521.21 


40. 


7.2± 0.7 


7.1 


6_2-7-i E 




85568.07 


75. 


5.6± 0.6 


4.7 


7 2 -6 3 A" 




86615.60 


103. 


3.3± 0.3 


3.0 


7 2 -6 3 A+ 




86902.95 


103. 


2.8± 0.3 


3.1 


8_4-9-3 E 




89505.78 


171. 


3.1± 0.4 


3.1 


l -2i E v t 


=1 


93196.66 


303. 


1.8± 0.2 


1.3 


8 3 -9 2 E 




94541.81 


131. 


4.8± 0.5 


5.0 


8 -7i A+ 




95169.52 


84. 


11. 5± 1.2 


8.6 


2i-li A+ v t 


= 1 


96396.06 


332. 


1.4± 0.2 


1.0 


2i-li E v t 


=1 


96492.16 


298. 


1.7± 0.2 


1.3 


2 -lo E v t - 


=1 


96493.55 


308. 


2.0± 0.2 


1.6 


2_i-l_i E v t = 


=1 


96501.70 


420. 


0.8± 0.2 


0.6 


2 -lo A+ v t 


=1 


96513.67 


431. 


0.7± 0.2 


0.7 


2i-li A" v t 


=1 


96588.59 


332. 


1.3± 0.2 


1.0 


2_i-l_i E 




96739.36 


12.5 


6.8± 0.7 


6.4 


2 -lo A+ 




96741.38 


7.0 


9.2± 0.9 


7.2 


2o-lo E 




96744.55 


20. 


8.1± 0.8 


7.1 


2i-li E 




96755.51 


28. 


6.8± 0.7 


6.1 


6i-5o E v t 


=1 


99730.96 


340. 


2.4± 0.4 


2.5 


13 2 -12 3 E 




100638.87 


234. 


4.5± 0.5 


5.6 


13-3-12-4 E 




104060.72 


274. 


2.3± 0.8 


3.7 


11-1-10-2 E 




104300.40 


159. 


6.8± 0.7 


7.7 


13_ 2 -13i E 




104336.64 


237. 


3.4± 1.1 


2.4 


IO4-H3 A" 




104354.86 


208. 


3.7± 0.8 


4.4 


IO4-H3 A+ 




104410.49 


208. 


3.4± 0.6 


4.4 


13i-12 2 A+ 




105063.76 


224. 


7.1± 0.7 


7.3 


14_2-14i E 




105576.39 


270. 


3.2± 0.3 


2.8 


3i-4 A+ 




107013.81 


28. 


7.5± 0.8 


9.9 


15_2-15i E 




107159.92 


305. 


3.6± 0.4 


3.3 


Oo-l-i E 




108893.93 


13.1 


6.9± 0.7 


7.3 


14 s -15 4 E 




109138.69 


380. 


4.3± 0.8 


2.9 


I6-2-I61 E 




109153.21 


342. 


3.6± 0.4 


3.5 


7 2 -8i A+ 




111289.60 


103. 


8.1± 0.8 


8.1 


5i-4 2 E 




216945.56 


56. 


28.0± 7.0 


37.8 


6i-7 2 A" v t 


=1 


217299.16 


374. 


15.5± 3.9 


24.2 


8 -7i E 




220078.49 


97. 


41.4±10.4 


47.6 
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a Integrated over the velocity range from —54 to —44 km s 
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Fig. 1. — Emission from lines with upper state energies less than 100 K for which data are available 
down to spatial frequencies of 3 kA or less. Lines are presented in order of increasing upper state 
energy. The map center is a(J2000) = 2 h 27 m 03?87, <5(J2000) = 61°52'24'.'6, the nominal position 
of W3(OH). Sutton et al. (2001) found the 3 mm continuum from the W3(OH) UCHii region to 
peak CK'2 ± O'/l west of this position. The TW object is 5'/94 east of the phase center (Wyrowski 
et al. 1999; Reid et al. 1995). With the exception of panel a, the plots are of average intensity over 
a v; sr range of —53 to —42 km s -1 . Panel a is a blend of lines, and the intensity is an average over 
a range of channels corresponding to —52.3 km s -1 for the 2o~lo A + line to —41.6 km s -1 for the 
2_i— l_i E line. Beams are shown in the lower left corners of the plots, and their sizes are listed in 
Table 1. Contours are .05, .1, .15, .2, .25, .3, .4, .5, .6 Jy beam -1 for panels a, c, d, e, and f and 
.075, .15, .225, .3, .4, .5, .6 Jy beam -1 for panel b. 

Fig. 2. — J = 2-1 lines from the first torsionally excited state. Map center and velocity range are 
the same as for Figure 1. Beams are shown in the lower left corners of the plots, and their sizes are 
listed in Table 1. Symbols mark the positions of the TW object and the 3 mm continuum center 
of W3(OH). Contours are multiples of 0.05 Jy beam -1 . 

Fig. 3. — Maps of 15 b-type lines towards W3(OH)/W3(H20). Map center and velocity range are 
the same as for Figure 1. Beams are shown in the lower left corners of the plots, and their sizes 
are listed in Table 1. Where space allows, symbols mark the positions of the TW object and the 
3 mm continuum center of W3(OH). Contours are .05, .1, .15, .4, .5, .6, .7, .8 Jy beam -1 for 
panels e, f, g, i, k, 1, m, and n. Contours are .125, .25, .375, .5, .75, 1., 1.5, 2. and then multiples 
of 1. Jy beam -1 for panels a, c, h, and j. For panels b, d, and o, the contours are multiples of 0.25 
Jy beam -1 . 

Fig. 4. — Continuum map at 1.3 mm with 0.6 arcsec resolution. Contours are multiples of 0.1 
Jy beam -1 . Peaks A, B, and C of Wyrowski et al. (1999) are labelled. Our measured position 
and positional uncertainty of the 3 mm W3(OH) continuum center is indicated by the filled box. 
The 2.0 x 1.2 arcsec ellipse represents our best estimate of the position and extent of the region of 
methanol emission near W3(H2 0). 

Fig. 5. — Rotation diagram fit to methanol emission from W3(H\20). In panel (a) the ordinate 
is 3c 2 / S v dv j (167r 3 z/ 3 f2// 2 £). In panel (b) the ordinate is the above quantity times C T , that is, 
corrected for optical depth. The solid line is the best weighted x 2 fit to the data, and gives a 
rotational temperature of 149 K and No/go = 1-3 x 10 15 cm -2 . In panel (c) the ordinate is the 
ratio of the observed integrated flux to that predicted by the non-LTE model discussed in the text. 

Fig. 6. — Spectra of three lines from a 4 x 4 arcsec region around the center of W3(H20), offset by 
(Aq = +5.3, A<5 = +0.1) arcsec from the reference position. 
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